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ABSTRACT 


This  thesis  addresses  the  optimization  of  all  types  of  space  electrical  propulsion 
thrusters.  From  the  Langmuir-Irving  payload  mass  fraction  formulation,  a  “dual¬ 
optimum”  solution  is  defined,  yielding  a  minimum  overall  mass  for  a  specified  payload 
consistent  with  minimum  transfer  time.  This  solution  fixes  the  ideal  payload  mass  ratio 
(mpi  /  mo)  at  a  value  of  0.45,  establishing  the  ratios  of  effective  exhaust  velocity  (v  /  vc) 
and  incremental  change  of  vehicle  velocity  (Au  /  vc)  to  characteristic  velocity  at  0.820 
and  0.327  respectively.  The  characteristic  velocity  (vc)  includes  thrust  time  as  well  as 
engine  efficiency  (rjt)  and  specific  power  (a).  A  range  of  mass  ratios  from  0.35  to  0.55  is 
used  in  order  to  allow  the  system  designer  some  flexibility  while  remaining  close  to 
optimal.  Nine  examples  are  presented  which  demonstrate  that  mission  profiles  can  be 
optimized  by  profile-to-thruster  matching.  A  comprehensive  list  of  currently  available 
electric  propulsion  engines  is  provided.  This  list  details  important  parameters  such  as  the 
specific  power,  which  “sizes”  an  engine  in  terms  of  power  provided  to  the  thruster  at  the 
cost  of  additional  mass.  Allowance  is  also  made  for  a  fuel  tank  mass  penalty,  and 
examples  show  that  this  can  also  noticeably  influence  the  optimum  design. 
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I.  INTRODUCTION 

A.  SPACE  ELECTRIC  PROPULSION  COMES  OF  AGE 

The  application  of  electric  propulsion  for  space  is  not  a  recent  revolution  in  the 
field.  The  idea  was  proposed  as  early  as  the  1950’s.  One  significant  article  appeared  as  a 
chapter  in  the  text,  Space  Technology,  entitled  “Low  Thrust  Flight:  Constant  Exhaust 
Velocity  in  Field-free  Space,”  written  by  D.  B.  Langmuir  (Langmuir,  1959,  pg.  9-01). 
Aside  from  the  evolution  of  technology,  little  has  changed  in  the  fundamentals  of  electric 
propulsion  since  that  time. 

In  the  last  35  years,  more  than  300  electric  propulsion  thrusters  have  flown  on 
over  100  spacecraft,  and  in  the  next  decade,  a  significant  increase  in  electrical  propulsion 
employment  is  expected  (Filliben,  1997,  pg.  5).  Nonetheless,  electrical  propulsion  is 
only  now  seeing  rapid  and  wider-scale  introduction  into  spacecraft  propulsion  programs. 

The  first  electrical  propulsion  engines  flew  in  the  mid-1960’s.  The  former  Soviet 
Union  used  a  pulsed  plasma  thruster  (PPT)  to  provide  attitude  control  to  their  Zond-2/1 
satellite,  which  launched  in  1964  (Martinez-Sanchez  and  Pollard,  1998,  pg.  692).  Then, 
in  1965,  the  United  States  followed  up  their  effort  using  a  resistojet  engine  for  orbit 
adjustment  of  the  Vela/2  satellite  (Martinez-Sanchez  and  Pollard,  1998,  pg.  689). 

If  both  of  those  missions  proved  successful,  why  has  additional  implementation  of 
electrical  propulsion  for  space  applications  been  so  seemingly  slow,  spanning  almost  40 
years?  For  this,  there  are  both  technical  and  doctrinal  reasons. 

First,  engineers  have  only  recently  been  able  to  package  sufficient  electrical 
power  and  power  conditioning  equipment  onboard  spacecraft  in  order  to  provide  for  the 
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demands  of  the  electrical  propulsion  system.  Unlike  chemical  propulsion,  electrical 
requires  energy  accumulation  and  energy  processing  functions  creating  a  power  supply 
penalty  at  the  expense  of  payload  mass  (Humble,  Henry,  and  Larson,  1995,  pg.  510). 
Likewise,  many  mission  planners  have  been  reluctant  to  replace  proven  chemical  rockets 
with  comparatively  untested  electrical  systems.  As  a  result,  electrical  systems  have  only 
gradually  become  accepted  and  only  for  less  demanding  tasks  --  such  as  north-south 
station  keeping  (NSSK)  and  orbit  insertion. 

Nonetheless,  electrical  propulsion  has  several  advantages  over  chemical  rockets 
for  near-zero  gravity  applications.  Moreover,  for  deep-space  missions,  which  require 
huge  incremental  changes  in  vehicle  velocity  (the  factor  Au),  electrical  propulsion  is  the 
only  realistic  option.  In  fact,  it  is  likely  that  the  future  of  deep  space  exploration  will 
depend  largely  upon  electric  propulsion  systems  as  chemical  systems  will  simply  be 
unable  to  attain  the  high  Au  required  for  these  missions.  Table  1  (Typical  Mission 
Velocity  Requirements)  presents  a  summary  of  various  mission  Au  requirements. 


Earth  to  LEOa 

7,600  m/s 

LEO  to  GEOb 

4,200 

LEO  to  Earth  escape 

3,200 

LEO  to  lunar  orbit  (7  days) 

3,900 

LEO  to  Mars  orbit  (0.7  yr) 

5,700 

LEO  to  Mars  orbit  (40  days) 

85,000 

LEO  to  Neptune  orbit  (29.9  yr) 

13,400 

LEO  to  Neptune  orbit  (50  yr) 

70,000 

LEO  to  solar  escape 

8,700 

LEO  to  1000  AUC  (50  yr) 

142,000 

LEO  to  Alpha-Centauri  (50  yr) 

30  x  106 

a  Low  Earth  Orbit  =  270  km 
b  Geosynchronous  Earth  Orbit  =  42,227  km  radius 
c  Astronomical  Units  (1  AU  =  149,558,000  km) 

Table  1:  Typical  Mission  Velocity  Requirements.  From  Hill,  Peterson,  page  508. 
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Note  however,  that  these  mission  velocity  requirements  are  calculated  for  chemical 
propulsion  and  not  electric  propulsion  systems.  Regardless,  Table  1  illustrates  the  typical 
magnitude  of  Au  required  for  various  missions  as  the  numbers  do  not  change 
significantly.  The  values  will  be  slightly  higher  for  electrical  systems. 

On  the  down  side,  electrical  propulsion  systems  cannot  overcome  high-gravity 
fields  due  to  their  low  acceleration  or  low  thrust.  Therefore,  the  initial  boost  to  attain 
escape  velocity  from  Earth  will  need  to  be  accomplished  via  chemical  rockets,  and  it  is 
evident  that  the  types  of  missions  slated  for  the  future  will  be  accomplished  via  a 
combination  of  both  electrical  and  chemical  propulsion.  There  are  specific  performance 
limitations  to  each.  (Hill  and  Peterson,  1992,  pg.  491) 

Chemical  rockets  can  be  described  as  “energy  limited”  in  that  the  fundamental 
chemical  behavior  of  the  propellant  defines  the  quantity  of  energy  (per  unit  mass  of 
propellant)  that  can  be  released  during  combustion.  High  propellant  energy  for 
spacecraft  is  possible  if  a  separate  energy  source  is  utilized  --  such  as  nuclear  or  solar 
energy.  In  that  manner,  the  rate  of  conversion  of  nuclear  or  solar  energy  to  electrical  and 
then  to  propellant  kinetic  energy  is  limited  by  the  mass  of  conversion  equipment  required. 
That  mass  can  be  a  comparatively  large  portion  of  the  total  vehicle  mass.  As  a  result,  the 
electrical  rocket  is  generally  said  to  be  “power  limited.”  (Hill  and  Peterson,  1992, 
pg.  491) 

B.  PHYSICS  BEHIND  ELECTRIC  PROPULSION 

Space  electrical  propulsion  systems  consist  of  three  basic  subsystems.  Those 
subsystems  include  the  power  processing  unit,  propellant  management  assembly,  and  the 
thruster  assembly  consisting  of  one  or  more  thrusters. 
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The  power  processing  unit  (PPU)  contains  electrical  interfaces  that  link  the 
spacecraft  power  source  and  the  propulsion  system.  The  PPU  sends  power  and  command 
signals  to  the  thruster(s)  and  propellant  flow  control  signals  to  the  propellant  subsystem. 
Additionally,  the  PPU  is  responsible  for  converting  various  power  sources  (solar, 
batteries,  fuel  cells,  nuclear,  etc.)  to  the  proper  voltage,  frequency,  pulse  rate,  and  current 
suitable  for  the  system  components. 

The  propellant  management  assembly  (PMA)  provides  the  means  of  storing, 
metering  and  delivering  the  propellant.  The  thrusters  convert  electrical  energy  to  kinetic 
energy  of  a  propellant  in  the  form  of  exhaust  or  thrust. 

Regardless  of  type,  the  common  feature  of  electric  propulsion  systems  is  the 
addition  of  energy  to  a  working  fluid  from  an  electrical  source.  Operating  in  a  steady  or 
pulsed  mode,  the  exhaust  is  accelerated  via  one  or  more  of  the  three  fundamental  types  of 
thrusters.  Those  types  of  thrusters  are  (1)  electrothermal,  (2)  electrostatic,  and  (3) 
electromagnetic.  (Martinez-Sanchez  and  Pollard,  1998,  pg.  688) 

In  reading  the  following  sections,  reference  Figure  1,  “Operating  Principles  of 
Various  Electrical  Propulsion  Thrusters,”  found  at  the  end  of  this  chapter  on  page  14. 

The  figure  will  prove  exceptionally  useful  in  understanding  the  operation  of  several  of 
the  following  thrusters  as  they  are  discussed. 

1.  Electrothermal  Thrusters 

In  the  electrothermal  thruster,  the  propellant  is  electrically  heated  and 
thermodynamically  expanded.  The  heated  gas  then  achieves  supersonic  velocity  as  it 
passes  through  a  nozzle.  Overall,  these  thrusters  have  one  major  drawback  in  that,  as  a 
class,  they  produce  limited  high  exhaust  velocities.  Additionally,  their  material 
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characteristics  limit  their  performance  to  values  that  are  generally  in  line  with  those 
exhibited  by  chemical  rockets.  At  present,  there  are  two  basic  types  of  electrothermal 
thrusters  --  resistojets  and  arcjets. 

The  resistojet  is  the  more  elementary  of  the  two  in  its  class.  In  essence,  it 
operates  by  using  an  electric  current  to  directly  heat  metal  components  (coiled  wire, 
tubing,  or  fins)  which  then  transfer  heat  to  the  gaseous  propellant  through  radiation 
and/or  convection.  As  a  result,  the  resistojet’s  performance  is  limited  by  the  capacity  of 
its  structural  components  to  withstand  high  temperatures,  and  resistojets  achieve  only  a 
modest  specific  impulse  of  300-310  seconds.  This  is  due  to  the  relatively  high  molecular 
mass  of  the  propellant  gas  in  addition  to  material  limitations  of  the  heating  wall,  which 
can  sustain  temperatures  up  to  approximately  2000  K.  (Martinez-Sanchez  and  Pollard, 
1998,  pg.  689) 

Even  with  its  comparatively  lower  value  of  specific  impulse,  the  resistojet’s 
superior  efficiency  contributes  to  far  higher  values  of  thrust/power  than  other  varieties  of 
electrical  propulsion.  Additionally,  these  engines  possess  the  lowest  overall  system  “dry 
mass”  since  they  do  not  require  a  power  processor  and  their  plumes  are  uncharged. 
(Sackheim  and  Byers,  1998,  pg.  672) 

Resistojets  have  seen  successful  employment  on  such  craft  as  Intelsat  V,  Satcom 
1-R,  GOMS,  Meteor  3-1,  Gstar-3,  and  Iridium  S/C  in  addition  to  older  satellites  and  test 
flights  (Martinez-Sanchez  and  Pollard,  1998,  pg.  689).  They  are  most  attractive  to  the 
mission  planner  for  low-to-modest  energy  applications.  In  particular,  they  are  most 
suitable  where  power  limits,  thrusting  times,  and  plume  impacts  are  mission  drivers. 
(Sackheim  and  Byers,  1998,  pg.  672) 
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The  other  form  of  the  electrothermal  thruster,  the  arcjet,  manages  to  overcome  the 
thermal  problems  encountered  by  the  resistojet,  and  it  enables  a  significant  increase 
(approximately  double)  in  specific  impulse  over  the  resistojet.  The  arcjet  accomplishes 
this  by  forming  an  attached  arc  within  the  nozzle  that  directly  heats  the  propellant  stream 
to  temperatures  much  greater  than  those  encountered  at  the  thruster  body.  Essentially, 
the  power  is  deposited  internally  in  the  form  of  the  electric  arc,  which  typically  operates 
between  a  concentric  upstream  rod  cathode  and  a  downstream  anode  that  serves  as  a 
supersonic  nozzle.  (Martinez-Sanchez  and  Pollard,  1998,  pg.  689) 

The  arc  core  reaches  temperatures  of  10-20,000  K  while  the  buffer  layer  near  the 
wall  is  maintained  at  less  than  2000  K.  As  a  consequence  of  the  temperature  gradient, 
there  is  practically  no  propellant  flow  through  the  arc  core.  This  results  in  high  specific 
impulse  at  the  cost  of  reduced  propulsive  efficiency.  In  addition,  the  flow  structure  at  the 
throat  is  extremely  non-uniform,  and  the  thruster  is  primarily  limited  by  arc  instabilities 
and  erosion  of  the  nozzle.  (Martinez-Sanchez  and  Pollard,  1998,  pg.  689)  Generally, 
arcjets  exhibit  significant  decreases  (approximately  six  times)  in  thrust/power  relative  to 
the  resistojet.  This  is  due  to  an  increased  specific  impulse  coupled  with  relatively  low 
values  of  efficiency  ranging  from  0.3  to  0.4.  (Sackheim  and  Byers,  1998,  pg.  672) 

One  significant  disadvantage  to  the  arcjet  is  that  the  required  PPU  is  substantially 
more  complex  than  that  encountered  in  the  resistojet  as  the  complex  plasma  arc  must  be 
controlled.  In  fact,  the  PPU  may  be  several  times  heavier  than  the  thruster  itself.  Still, 
however,  the  arcjets  exhibit  relatively  low  dry  masses  when  compared  with  some  of  the 
other  electrical  propulsion  systems  in  use  or  development.  (Martinez-Sanchez  and 


Pollard,  1998,  pg.  689) 


With  its  intermediate  specific  impulse  of  600  -  650  seconds,  the  arcjet  is  a  very 
viable  option  for  short  bum  duration  missions.  The  lower  specific  impulse  implies  higher 
thrust  for  a  given  power,  and  arcjets  have  been  used  for  geostationary  applications, 
including  the  Telstar  4  and  GE-1  satellites.  (Martinez-Sanchez  and  Pollard,  1998,  pg. 
689)  The  engines  are  relatively  simple  to  integrate,  and  they  are  the  least  costly  and 
complex  of  any  plasma  propulsion  device  (Sackheim  and  Byers,  1998,  pg.  672). 

2.  Electrostatic  Thrusters 

Electrostatic  thrusters  rely  on  electrostatic  fields  in  order  to  exert  forces  on  the 
propellant  which  is  throttled  into  the  field  in  the  form  of  charged  particles.  Such  thrusters 
can  only  operate  in  a  vacuum,  and  all  the  created  particles  must  be  of  the  same  sign  —  i.e. 
all  positive  or  all  negative.  Likewise,  those  particles  must  be  neutralized  after  passing 
through  the  engine.  Otherwise,  the  thruster  would  eventually  render  itself  ineffective  as  a 
net  charge  would  build  up  on  the  spacecraft  setting  up  the  return  of  oppositely  charged 
ions  to  the  craft.  As  a  result,  any  thrust  would  be  cancelled.  In  addition,  sensitive  system 
components  could  be  damaged  by  the  returning  ions.  (Humble,  Henry,  and  Larson,  1995, 
pg.  532) 

Naturally,  electrons  would  seem  to  be  the  most  likely  candidate  particles  for  these 
engines.  They  are  easy  to  produce  and  to  accelerate;  however,  their  small  mass  renders 
them  impractical.  They  are,  in  fact,  so  extremely  light  that  the  momentum  imparted  to 
them  is  negligible  even  with  their  high  velocity. 

Consequently,  electrostatic  thrusters  operate  using  heavy  atoms  charged  as 
positive  ions  or  as  colloids  (liquid  droplets).  Neutralization  of  the  positive  ions  is  easily 
accomplished  with  electrons,  and  the  ions  can  be  240,000  times  as  heavy  as  the  electron. 
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The  charged  colloids  are  generally  at  least  10,000  times  as  massive  as  the  ions.  Thus, 
substantial  momentum  can  be  generated. 

The  implementation  of  heavy  atoms  has  other  advantages.  The  thrust  per  unit 
area  increases  as  the  square  of  the  particle  mass  to  the  charge  ratio  in  addition  to  the 
highly  desirable  characteristic  of  high  voltage  and  low  current  as  opposed  to  low  voltage 
and  high  current  (Sutton,  1992,  pg.  580). 

Electrostatic  thrusters  are  categorized  by  the  source  of  charged  particles,  and  all 
three  types  exhibit  the  same  three  basic  stages:  ion  production,  acceleration,  and 
neutralization.  The  five  types  of  thrusters  that  have  been  studied  are  as  follows: 

1 .  Electron  Bombardment  Ion  Thruster.  Positive  ions  are  produced  via  the 
bombardment  of  vaporized  or  gaseous  propellant  (such  as  xenon  or  mercury) 
with  electrons  emitted  from  a  heated  cathode. 

2.  Ion  Contact  Thruster.  Positive  ions  are  created  by  passing  a  propellant  vapor 
(cesium)  through  a  hot  contact  ionizer  (usually  tungsten).  Although  it  was 
intensely  investigated  15  to  30  years  ago,  the  ion  contact  thruster  was 
abandoned  as  impractical. 

3.  Field  Emission  (Colloid)  Thruster.  This  engine  generates  positive  or 
negative  particles  by  passing  tiny  liquid  droplets  of  propellant  through  an 
intense  electric  field  (corona  discharge).  (Sutton,  1992,  pg.  580) 

4.  Radio  Frequency  Ion  Thruster.  This  thruster  produces  a  high  specific  impulse 
at  the  price  of  complexity  and  mass.  Essentially,  it  consists  of  a  discharge 
chamber  that  is  surrounded  by  an  inductive  coil  connected  to  a  radio 
frequency  generator.  Neutral  propellant  atoms  are  ionized  by  the 
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bombardment  of  electrons  which  are  accelerated  by  the  induced  high- 
frequency  eddy  fields  in  the  discharge  chamber. 

5.  Electron  Cyclotron  Resonance  /  Microwave  Discharge  Ion  Thrusters.  Here  a 
circularly  polarized  microwave  beam  is  projected  from  a  conventional 
waveguide  via  a  suitable  dielectric  window  in  order  to  ionize  the  propellant  in 
the  discharge  chamber.  (Filliben,  1997,  pg.  5) 

At  the  present  time,  ion  engines  are  the  “thruster  of  choice”  for  deep  missions 
such  as  interplanetary  transfers.  These  missions  require  a  high  Au  and  must  tolerate  long 
thrusting  times.  (Martinez-Sanchez  and  Pollard,  1998,  pg.  692)  A  typical  xenon  ion 
thruster  can  generate  an  average  mission  specific  impulse  of  2,800  -  3,500  seconds  with  a 
lifetime  of  greater  than  8000  hours  (333  days). 

3.  Electromagnetic  Thrusters 

There  are  several  variations  to  the  electromagnetic  type  of  thruster.  All  make  use 
of  a  propellant  gas  that  is  heated  to  a  plasma  state.  That  plasma  is  then  part  of  a  current- 
carrying  electric  circuit  which  interacts  with  combined  electric  and  magnetic  fields 
thereby  generating  thrust. 

The  four  types  of  electromagnetic  thrusters  which  we  will  discuss  include  the 
Hall-effect  thruster,  stationary  plasma  thruster  (SPT),  the  magnetoplasmadynamic  (MPD) 
thruster,  and  the  pulsed  plasma  thruster  (PPT). 

The  Hall-effect  thruster  is  named  for  its  closed  circular  electron  drift  that  is 
exhibited  between  the  cathode  and  anode  of  the  thruster.  This  effect  is  analogous  to  the 
Hall  Effect.  The  SPT  operates  on  similar  principles. 
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In  the  Hall-effect  thruster,  an  axial  electric  field  and  radial  magnetic  field  are 
established  in  the  discharge  chamber.  Meanwhile,  neutral  propellant  atoms  are  fed  into 
the  discharge  chamber  and  ionized  by  the  bombardment  of  electrons.  The  perpendicular 
electric  and  magnetic  fields  then  force  the  electrons,  emitted  from  an  external  hollow 
cathode,  to  become  trapped  in  azimuthal  drift  motion.  At  the  same  time,  the  radial 
magnetic  field  is  not  sufficiently  strong  to  affect  the  propellant  ion  trajectory,  and  they 
are  consequently  accelerated  axially  by  the  electric  field,  producing  thrust.  (Filliben, 
1997,  pg.  5) 

The  first  use  of  a  Hall  thruster  electric  propulsion  system  in  space  took  place  in 
October  of  1998.  The  U.  S.  Navy’s  Hall  thruster,  installed  on  an  experimental  spacecraft 
owned  by  the  National  Reconnaissance  Office,  began  operating  as  part  of  the  electric 
propulsion  demonstration  module  (EPDM).  The  engine  provides  for  orbit  raising  and 
long-term  station  keeping.  (Space  Tracks,  1999,  pg.  13) 

The  particular  engine  onboard  the  EPDM  is  a  TAL  (thruster  with  anode  layer) 
type  Hall  thruster,  using  xenon  gas.  It  is  a  gridless  device  producing  greater  than  twice 
the  thrust  of  an  ion  engine  of  the  same  operating  power.  (Space  Tracks,  1999,  pg.  13) 
Although  still  in  development,  the  magnetoplasmadynamic  thruster  (MPD)  is 
generally  regarded  as  the  leading  potential  candidate  for  future  deep  space  missions  such 
as  a  heavy-lift  Mars  transfer,  utilizing  a  nuclear  powerplant.  (Martinez-Sanchez  and 
Pollard,  1998,  pg.  693) 

The  engine  operates  by  generating  a  current  along  a  conducting  bar.  The  current 
consequently  yields  a  self-induced  azimuthal  magnetic  field  that  interacts  with  the  current 
of  an  arc  that  travels  from  the  point  of  the  bar  to  a  conducting  wall.  The  resulting  force 
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has  two  components.  One  component  is  a  radially  inward  force  that  constricts  the  flow, 
and  the  other  force  acts  along  the  axis,  producing  directed  thrust.  As  expected,  one  of  the 
great  problems  with  this  thruster  is  the  substantial  erosion  that  occurs  at  the  point  of 
contact  between  the  arc  current  and  electrodes.  (Santarius,  1997)  Additionally,  the  MPD 
thruster  is  more  difficult  to  optimize  due  to  the  fact  that  several  physical  effects  are 
involved  in  its  operation. 

The  Pulsed  Plasma  Thruster  (PPT)  is  “different  from  all  other  concepts”  in  two 
fundamental  ways.  This  device  operates  in  short  pulses  of  approximately  10 
microseconds  duration,  and  in  its  most  developed  form  uses  a  portion  of  solid  propellant 
feedstock  (typically  Teflon®).  As  a  result,  no  propellant  tanks  (referred  to  as  “tankage”) 
are  required  in  this  device.  (Martinez-Sanchez  and  Pollard,  1998,  pg.  692) 

In  the  PPT,  a  portion  of  the  feedstock  is  ablated  and  ionized  by  an  electrical  arc 
discharge  sheet  initiated  between  two  electrodes  by  a  discharging  capacitor.  The 
resulting  propellant  plasma  is  accelerated  by  the  interaction  of  the  arc  and  the  self- 
induced  magnetic  field  of  the  current  loop.  (Filliben,  1997,  pg.  5) 

The  PPT  design  is  comparatively  simple  and  reliable.  However,  one  of  the 
inherent  problems  with  this  thruster  is  that  the  force  falls  off  as  the  current  loops  get 
large.  In  addition,  PPT’s  generally  suffer  from  low  efficiencies  on  the  order  of  8  to  13 
percent.  (Santarius,  1997) 

PPT’s  offer  two  additional  distinct  advantages  in  that  they  integrate  a  non-toxic 
propellant  feed  system  with  a  thruster  in  a  single  compact  unit  as  well  as  offering  a 
variable  pulse  repetition  rate.  As  a  result,  they  possess  a  flexibility  of  operation  over  a 
wide  range  of  mean  power  or  thrust,  and  therefore,  they  have  found  employment  for 
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precision  orbital  or  attitude-control  tasks  since  the  late  1960’s  with  the  LES-6  satellite 
and  the  U.  S.  Navy’s  NOVA  constellation.  (Martinez- Sanchez  and  Pollard,  1998, 
pg.  692) 

The  optimal  implementation  of  these  three  fundamental  types  of  electric 
propulsion  engines  (electrothermal,  electrostatic,  and  electromagnetic)  will  be  the  focus 
of  our  analysis.  As  illustrated  in  this  chapter,  there  are  several  variants  of  each  of  the 
engines.  Some  of  those  engines  are  currently  in  use  or  in  various  stages  of  development. 
Please  reference  Table  2  on  page  13  for  a  summary  of  typical  electrical  propulsion 
performance  parameters. 
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Type  Thrust  (mN)  /  Duration  Is(sec)  Efficiency3  Specific  Power  (W/mN)  Typical  Propellants 
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II.  THEORETICAL  BACKGROUND 


A.  FIRST  PRINCIPLES 

Electrical  propulsion  systems  have  four  general  yet  unique  characteristics.  These 
characteristics  are  exploited  in  analyzing  and  understanding  the  electric  propulsion 
system,  and  they  are  delineated  in  the  following  paragraph. 

Electric  propulsion  engines  (1)  incorporate  a  comparatively  heavy  mass  of 
required  power  generation  equipment,  and  they  (2)  produce  low  thrust.  As  a  result  of  that 
low  thrust,  (3)  very  low  accelerations  (on  the  order  of  lCT4  to  10'6  g0)  are  attained  (4) 
requiring  long  operating  times.  In  the  following  development  of  the  equations  governing 
electrical  engines,  the  electrical  rocket’s  performance  will  be  characterized  in  terms  of 
power  and  mass  as  first  proposed  by  Langmuir  and  Irving.  (Sutton,  1992,  pg.  596) 

First,  the  mass  aspects  of  electrical  propulsion  will  be  analyzed. 

The  major  components  integral  to  an  electrically  powered  space  vehicle  are 
broken  down  as  follows: 

mo  =  mp  +mpi  +mpp.  (1) 

In  words,  this  equation  states  that  the  initial  total  mass  of  the  vehicle  (m<j)  is  equivalent  to 
the  sum  of  the  masses  of  the  propellant  (mp),  payload  (mpi),  and  powerplant  (mpp).  Refer 
to  the  list  of  symbols  and  abbreviations  (on  page  xiii)  for  the  comprehensive  table  of  the 
following  abbreviations  and  their  respective  units. 

The  determination  of  mpp  for  an  assigned  mission  is  the  most  amorphous.  This 
term  incorporates  a  wide  variety  of  components  that  make  up  the  empty  (or  un-foeled) 
propulsion  system.  Those  components  include  the  thruster,  the  propulsion  storage  and 
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feed  system,  the  energy  source  with  its  conversion  system  and  auxiliaries  as  well  as  the 
associated  structure  (Sutton,  1992,  pg.  597). 

Now,  let  us  look  at  the  power  aspects. 

Regardless  of  the  energy  source  (battery,  nuclear,  solar,  fuel  cells,  etc.)  to  the 
input  power  supply,  the  input  energy  must  always  be  greater  than  the  required  electrical 
power  output,  and  power  conversion  efficiencies  typically  average  approximately  70 
percent.  This  aspect  is  accounted  for  in  the  engine  parameter  known  as  specific  power, 
a,  or  the  “state-of-the-art.”  It  is  defined  as  the  ratio  of  electrical  power  output  to  the  mass 
of  the  power  plant,  and  it  serves  as  an  indication  of  the  system’s  ability  to  effectively  and 
efficiently  generate  or  convert  to  electrical  power.  The  equation  for  a  is  defined  as 
follows: 

a  =  Pe/ mpp,  (2) 

where  Pe  is  the  electrical  power  that  is  actually  supplied  to  the  propulsion  system. 

Typical  high  “alpha”  values  that  are  currently  attainable  range  from  approximately  100  to 
200  W/kg.  It  is  a  singular  goal  in  electrical  propulsion  to  strive  to  push  this  number  to 
the  highest  attainable  values.  (Sutton,  1992,  pg.  597) 

In  an  electrical  propulsion  thruster,  the  electrical  power  input  is  converted  to  the 
kinetic  energy  of  the  propellant  exhaust.  Thus,  Pe  can  also  be  defined  by  accounting  for 
losses  due  to  thruster  efficiency,  qt,  where 

Tit  =  Pjet  /  Pe  =  (1/2)  m  V2  /  Pe  (Sutton,  1992,  pg.  571),  (3) 

and,  consequently,  we  may  redefine  Pe  as  the  following: 

Pe  =  a  mpp  =  (1/2)  m  v2  /  r|t  =  mp  v2  /  (2  tp  Tit).  (4) 
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In  this  equation,  m  is  the  propellant  mass  flow  rate,  v  is  the  effective  exhaust  velocity, 
and  tp  is  the  duration  of  operation  in  which  the  propellant  mass,  mp,  is  ejected  from  the 
thruster.  It  is  important  to  note  that  this  equation  assumes  a  uniform  rate  of  propellant 
expulsion.  (Sutton,  1992,  pg.  597) 

Using  the  previous  definitions  of  nio,  a,  and  Pe,  we  can  form  the  equation  for  the 
payload  mass  fraction,  which  is  written, 

mpi  /  m0  =  (1  -  (v2/(2  a  tp  r|t))(eAu/v  - 1))/  eAu/v.  (5) 

This  equation  assumes  gravity-free  and  drag-free  flight.  Additionally,  by  defining  the 
characteristic  velocity,  vc,  as 

vc2  =  2  a  tp  r|t,  (6) 

we  can  re-write  the  payload  mass  fraction  as  follows: 

mp,  /  mo  =  (1  -  (v  /  vc)2(eAu/v  - 1))/  eAu/v.  (7) 

In  this  form,  the  payload  mass  fraction  is  more  manageable.  (Sutton,  1992,  pg.  597) 

A  plot  of  this  equation  is  illustrated  in  Figure  2.  The  equation  is  graphed  over  the 
range  of  v  /  vc  for  various  values  of  Au  /  ve.  It  is  evident  that  a  maximum  exists  such  that 
the  payload  mass  ratio,  mpi  /  mo,  for  a  given  velocity  increment,  Au  /  vc,  is  maximum  at  a 
particular  value  of  exhaust  velocity,  v  /  vc,  or  specific  impulse,  Is  /  vc.  Specific  impulse 
and  effective  exhaust  velocity  are  related  via  the  equation, 

Is =  v  /  g0  where  g0=  9.81  m  /  s2.  (Sutton,  1992,  pg.  597)  (8) 

The  maximum  described  in  these  curves  exists  for  two  reasons.  First,  the  inert 
mass  of  the  power  plant,  mpp,  increases  almost  linearly  with  the  specific  impulse. 
However,  the  mass  of  consumed  propellant  decreases  with  specific  impulse.  The  vehicle 
velocity  increment,  Au,  increases  with  higher  Is;  however,  Au  is  also  decreased  by  the 
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lower  mass  ratio  due  to  the  increased  inert  mass.  Consequently,  there  is  a  unique 
optimum  value  of  Is  (or  v)  for  every  electric  propulsion  mission  as  defined  by  the 
corresponding  required  value  of  Au.  (Sutton,  1992,  pg.  597) 

Again,  refer  to  Table  1  for  a  brief  summary  of  various  mission  Au  values. 


Figure  2.  Payload  Mass  Fraction  Plot:  mpi  /  nio  vs.  v  /  vc. 

The  payload  mass  ratio  can  also  be  re-written  to  a  somewhat  more  useful  form  by 
solving  for  Au  /  vc.  In  this  case,  the  equation  takes  on  the  following  form: 

Au  /  vc  =  (v  /  vc)  In  [(1  +  (v  /  vc)2)/((mpi  /  mo)  +  (v  /  vc)2)],  (9) 

and  it  is  plotted  in  Figure  3.  In  this  case,  the  equation  is  graphed  over  the  same  range  of 
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A/ny 


v  /  vc  for  various  values  of  mpi  /  mo. 


This  graph  shows  the  optimum  curves  for  three  different  mass  payload  ratios. 
From  this  plot  and  the  defining  equations,  we  note  that  the  optimal  values  for  Au  /  vc  and 
v  /  vc  are  as  follows: 


mpi/nio 


Au/  vc 


0.35  0.404  0.782 

0.45  0.327  0.820 

0.55  0.257  0.864 


Table  3.  Optimum  Values:  Au  /  vc  and  v  /  ve  for  Various  Payload  Mass  Ratios. 
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These  values  will  be  an  essential  aspect  of  our  engine  selection  (or  engine-mission 
matching)  algorithm,  which  is  developed  in  Chapter  IV. 

B.  ANALYSIS  -  DUAL  OPTIMUM  METHODOLOGY 

Incorporating  the  previous  development  via  “first  principles,”  it  is  now 
advantageous  to  establish  a  “joint”  or  “dual”  optimum,  which  will  optimizes  the  payload 
mass  ratio  and  yields  the  minimum  transfer  time.  It  can  be  calculated  by  multiplying  the 
payload  mass  ratio  and  the  previously  defined  equation  for  Au  /  vc. 

The  plot  of  the  resulting  equation  (for  various  payload  ratios)  is  illustrated  in 
Figure  4,  and  it  shows  that  the  dual-optimum  payload  mass  ratio  peaks  at  a  value  of 
mpi  /  m0  =  0.45. 


Figure  4.  Dual  Optimum. 


20 


This  plot  indicates  that  the  optimum  payload  ratio  is  0.45.  However,  due  to  the 
fact  that  it  may  be  difficult  for  the  designer  to  achieve  this  mass  ratio,  an  approximate 
design  range  of  payload  mass  ratios  of  0.35  to  0.55  is  investigated.  As  illustrated  in 
Figure  4,  these  values  do  not  lie  too  far  below  the  optimum  for  0.45. 

It  should  be  restated  that  for  constant  propellant  mass  flow  (m)  constant 
thrust  (F)  and  negligibly  short  start  and  stop  transients: 

Is  =  F /  m  go  and  .\F  =  m  g0 Is.  (10) 

Hence,  for  a  higher  specific  impulse,  Is,  less  fuel  is  required.  However,  as  the  specific 
impulse  is  increased,  we  may  drive  through  the  optimum  due  to  the  equation  for  Au  /  vc, 
which  incorporates  the  applicable  “efficiencies”  of  the  engines  in  the  forms  of  a,  tp,  and 
rjt-  As  a  result,  an  engine  with  an  arbitrarily  high  Is  may  not  be  optimum  for  a  particular 
application  or  mission.  For  example,  the  Teflon  PPT  has  relatively  high  values  for  Is  at 
the  expense  of  comparatively  poor  values  for  a  and  r|t. 

Figure  5  illustrates  a  three-dimensional  optimum  surface.  The  plot  is  illustrated 
across  the  axis  values  of. . . 


x:  v/vc 

range: 

0.5  to  1.3 

y:  mpi/mo 

range: 

0.3  to  0.6 

z:  (mpi  /  nio)(Au  /  vc) 

range: 

0.1 15  to  maximum  at  0.148 

The  plot  clearly  shows  the  establishment  of  a  optimum  or  maximum  point  at  the  payload 
mass  ratio  of  precisely  mpi  /  mo  =  0.45. 
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Figure  5.  Three-dimensional  Dual  Optimum  Surface. 


C.  TANKAGE  PENALTY 

The  parametric  curves  that  have  been  presented  up  to  this  point  are  inherently 
ideal  in  nature.  They  generally  assume  that  all  spacecraft  mass  can  fall  neatly  into  the 
three  categories  of  mass  of  propellant,  payload,  or  power-plant.  However,  as  previously 
noted,  it  can  be  difficult  to  assess  where  exactly  certain  masses  (such  as  shared  power 
conditioning  equipment)  should  be  charged,  and  as  a  result,  the  calculation  of  the  specific 
power,  a,  for  the  respective  electric  propulsion  engine  can  have  high  variability.  In  the 
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same  manner,  one  can  assess  somewhat  of  a  “tankage  penalty”  in  proportion  to  the 
amount  of  propellant  required  as  a  result  of  the  following  reasoning. 

The  typical  propellant  tanks  are  composed  of  overwrapped  titanium  or  aluminum. 
Their  loading  varies  depending  upon  several  factors  including  satellite  mass,  expected 
life  on  orbit,  mission  (orbit  raising,  NSSK,  etc.)  and  cant  angle  of  thruster.  As  a  result, 
the  actual  size  and  mass  of  tanks  will  vary  regardless  of  engine  type  or  manufacturer. 

An  accurate  figure  of  merit  for  the  required  tankage  revolves  around  the  following 
assumptions.  First,  that  a  propellant  tank  that  holds  zero  propellant  mass  is  assessed  at 
one  kilogram.  Second,  that  the  tank  has  a  mass  fraction  of  approximately  ten  percent  as  a 
function  of  loaded  propellant.  It  is  important  to  note  that  all  of  the  listed  electrical 
propulsion  engines,  except  the  pulsed  MPD  type,  require  tankage.  (Clauss,  1999) 

For  example,  a  tank  that  holds  40  kg  of  xenon  would  have  a  dry  mass  of... 

1.0  kg  +  0.1  (40  kg )  =  5  kg. 

If  this  10  %  tankage  “penalty”  and  one  kilogram  of  dead  dry  mass  are 
incorporated  into  the  previous  analysis,  the  equation  for  the  initial  mass  and  components 
is  altered  from  mo  =  mp  +  mpi  +  mpp  to  mo  =  1 . 1  mp  +  mpi  +  mpp .  However,  the 
equation  for  final  mass  remains  unaltered  as  follows:  mf  =  mo  -  mp .  The  0.1  mp  is 
residual  mass  that  remains  within  the  system. 

Using  these  definitions,  the  equation  for  the  rewritten  payload  mass  ratio,  solved 
for  Au  /  vc,  takes  the  following  form: 

Au  /  vc  =  (v  /  ve)  In  [(1 . 1  +  (v  /  vc)2)/((mpi  /  nio)  +  (v  /  vc)2  +  0.1)].  (11) 

This  equation  is  plotted  as  Au  /  vc  vs.  v  /  vc  in  Figure  6. 
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Figure  6.  Payload  Mass  Fraction  Plot  with  10  %  Tankage  Penalty: 

Au  /  vc  vs.  v  /  vc. 

Note  that  the  curve  is  essentially  similar  to  the  plot  illustrated  in  Figure  3; 
however,  it  is  adjusted  for  the  ten  percent  tankage  “penalty.”  As  a  result,  the  optimum 
curves  have  shifted  down  and  right  from  their  original  values.  This  is  illustrated  in 
Figure  7  which  shows  the  optimum  curves  with  and  without  the  tankage  “penalty”  for  the 
payload  mass  ratio  (mpi  /  mo)  of  0.45. 

For  a  ten  percent  tankage  “penalty,”  the  following  optimal  values  for  Au  /  vc  and 
v  /  vc  result: 
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moi/mo 

Au  /  vc  |  opt 

v/vc  opt 

0.35 

0.375 

0.850 

0.45 

0.306 

0.895 

0.55 

0.240 

0.945 

Table  4.  Optimal  Values  with  10%  Tankage  Penalty:  Au  /  vc  and  v  /  vc. 


Figure  7.  Optimum  Comparison  with  and  without  10  %  Tankage  Penalty: 

Au  /  v0  vs.  v  /  vc  at  mpi  /  mo  =  0.45. 

For  all  of  the  plotted  “tankage”  penalty  curves,  the  values  plotted  are  for  a  10  % 
penalty.  Typical  tankage  fractions  range  from  0.05  for  the  hydrazine  (N2H4)  resistojet 
and  arcjet  to  0.1 5  for  the  H2  arcjet.  Xenon  Hall  thrusters  and  ion  thrusters  have  tankage 
fractions  of  about  0.12.  (Martinez-Sanchez  and  Pollard,  1998,  pg.  690) 
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III.  CURRENT  TECHNOLOGY  -  ELECTRIC  ENGINES 

A.  NASA  DEVELOPMENT  AND  OBJECTIVES 

Today,  there  are  many  organizations  pursuing  the  research  and  development  of 
space  electric  propulsion.  Various  private  companies  (including  Primex,  Hughes,  ARC, 
and  Daimler  Chrysler)  presently  offer  electric  propulsion  engines  for  sale  and  use  in 
various  applications.  Additionally,  several  research  facilities  —  such  as  NRL  (Naval 
Research  Labs)  and  NASA’s  JPL  (Jet  Propulsion  Labs)  --  have  also  developed  various 
engines  or  components  of  engines.  However,  the  one  department  of  NASA  that  is 
primarily  responsible  for  electric  propulsion  research  is  the  On-board  Propulsion  Branch. 

NASA’s  On-board  Propulsion  Branch  is  based  at  the  John  H.  Glenn  Research 
Center  located  on  Lewis  Field  in  Cleveland,  Ohio.  The  branch  is  responsible  for 
developing  advanced  on-board  spacecraft  propulsion  technologies  --  both  chemical  and 
electrical  —  for  future  NASA  missions.  The  electrical  propulsion  technologies  currently 
under  investigation  at  Lewis  Field  include  resistojets,  arcjets,  gridded  ion  xenon  thrusters, 
Hall  xenon  ion  thrusters,  and  pulsed  plasma  thrusters  (PPT’s).  (Oleson,  1999) 

In  October  of  1998,  the  on-board  propulsion  branch  celebrated  two  major 
milestones  in  the  realm  of  electrical  propulsion.  Deep  Space  One  (DS1)  was  successfully 
launched  on  October  24,  1998  from  Florida’s  Cape  Canaveral  Air  Station.  This  craft  is 
intended  to  validate  new  space  technologies,  including  its  ion  propulsion  primary  engine, 
which  was  developed  in  the  NSTAR  (NASA  Solar  Electric  Propulsion  Technology 
Application  Readiness)  program.  NASA  Glenn  developed  the  engine  and  power 
processors  for  DS1,  and  this  flight  was  the  first-ever  application  of  such  a  technology  to  a 
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deep  space  mission.  DS1  was  to  rendezvous  with  the  near-Earth  asteroid,  1992  KD,  in 
July  1999,  and  it  passed  within  15  miles  of  the  asteroid  Braille  on  July  28  at  a  speed  of 
greater  than  35,000  mph  (15.5  km/s).  (Oleson,  1999) 

DS 1  ’s  NSTAR  ion  engine  utilizes  81  kilograms  of  xenon  propellant,  which  it 
ionizes  and  accelerates  to  68,000  mph.  The  engine,  which  is  12  inches  in  diameter,  can 
operate  at  half  throttle  for  greater  than  20  months,  generating  a  continuous  thrust  of 
0.02  lbf.  As  a  result,  the  engine  is  ten-times  more  fuel  efficient  than  a  comparable 
chemical  system,  and  it  is  the  prototype  for  smaller  /  lower-cost  space  vehicles.  (Oleson, 
1999) 

On  October  3,  1998,  a  STEX  spacecraft  was  launched,  carrying  the  EPDM 
(Electric  Propulsion  Demonstration  Module).  This  module  was  developed  under  joint 
effort  on  the  part  of  NASA,  Naval  Research  Labs  (NRL)  and  private  industry.  The  Glenn 
Research  Center  led  the  program,  which  was  intended  to  demonstrate  the  first  use  of  a 
high  specific  impulse,  low-power  Hall  thruster  system  on  a  U.  S.  spacecraft.  (Oleson, 
1999) 

That  particular  Hall  engine  is  a  TAL  (Thruster  with  Anode  Layer)  D-55  Hall- 
effect  thruster.  It  was  developed  by  the  Russian  Central  Research  Institute  of  Machine 
Building  (TSNIIMASH),  and  it  successfully  operated,  firing  ten  times  from  23  to  24 
October.  The  engine  provided  40  mN  of  thrust  and  raised  the  satellite’s  orbit  0.35  NM. 
(Oleson,  1999) 

As  electrical  propulsion  evolves,  NASA  Glenn  has  clearly  defined  its  objectives. 

For  fiscal  year  1999,  four  major  product  targets  have  been  defined.  Those  targets 
include  a  next  generation  ion  system  for  near-  and  further-term  space  science  and  a 
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pulsed  plasma  thruster  (PPT)  technology  for  precision  imaging.  In  addition,  NASA 
would  like  to  develop  both  a  high  specific  impulse,  low-mass  Hall  system  and  a 
miniature  electrothermal  thruster  for  Earth  science.  (Oleson,  1999) 

NASA  Glenn  has  also  defined  what  it  calls  “grand  challenges”  in  the  development 
of  electrical  propulsion.  Those  challenges  are  delineated  for  each  of  the  three  types  of 
electrical  propulsion  systems,  and  they  are  outlined  as  follows:  (Oleson,  1999) 

Electrothermal  Thrusters: 

•  integrated  electronics  /  thruster  package  for  micro-spacecraft  applications 

Electrostatic  Thrusters: 

•  lightweight,  high  total  impulse  ion  systems  for  space  science 

•  very  high  total  impulse  ion  system  for  far-term  space  science 

Electromagnetic  Thrusters: 

•  low-cost,  simple,  low-power  Hall  for  miniature  Earth-orbital  spacecraft 

•  high  specific  impulse  (Is),  direct  drive  Hall  for  space  science  and  Earth  orbital 

•  increase  thruster  efficiency  by  factor  of  two  to  16% 

•  increase  capacitor  life  to  40,000,000  pulses  and  beyond 

•  lower  electromagnetic  emissions  to  below  present  standard  specifications 

•  predict  contamination  on  spacecraft  surfaces  in  design  phase 

Nonetheless,  the  ability  and  timeliness  of  electric  propulsion  to  meet  these 

challenges  is  unpredictable.  This  is  due  inherently  to  the  fact  that  the  “state  of  the  art”  in 
electric  propulsion  is  due  and  determined  largely  by  the  “state  of  the  art”  of  the  electrical 
components  that  comprise  various  aspects  of  the  electrical  propulsion  system.  Primarily, 
the  concern  is  the  required  additional  mass  of  electrical  components.  This  fact  will 
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become  evident  in  the  forthcoming  discussion  of  the  specific  power  or  “alpha”  (a) 
parameter. 

B.  SUMMARY  OF  AVAILABLE  ENGINE  DATA 

One  of  the  contributions  of  this  thesis  is  to  obtain  and  summarize  the  most  current 
reported  data  regarding  present-day  electric  propulsion  engine  performance 
characteristics.  This  summary,  which  includes  data  from  various  vendors  as  well  as  the 
above  research  sites,  can  be  found  on  page  36  in  Table  5:  Summary  of  Current 
Technology  Engines,  and  the  various  sources  of  this  information  are  listed  at  the  bottom 
of  the  table.  Figure  8,  on  page  35,  provides  an  overall  map  of  current  engines  and  their 
respective  values  of  specific  impulse,  Is. 

The  electrical  engine  parameter  that  has  certainly  been  the  most  difficult  to  define 
has  been  an  accurate  representation  of  the  parameter  a,  specific  power  or  “state  of  the 
art.”  This  parameter  is  so  named  because  it  defines  the  relation  of  the  power  conversion 
unit  with  respect  to  the  system’s  mass.  In  other  words,  as  previously  defined, 
a  =  Pe  /  mpp ,  where  Pe  is  the  power  supplied  to  the  electrical  propulsion  system  and  mpp 
is  the  power-plant  mass  (Sutton,  1992,  pg.  597). 

The  apparent  controversy  surrounding  the  “state  of  the  art”  parameter  is 
evidenced  in  the  article  “Optimization  of  Electric  Propulsion  Systems  Considering 
Specific  Power  as  Function  of  Specific  Impulse”  by  M.  Auweter-Kurtz  et  al.  This  article, 
found  in  the  Nov.-Dee.  1988  issue  of  the  AIAA  Journal  of  Propulsion  and  Power, 
presents  the  argument  that  specific  power,  a,  is  a  strong  function  of  specific  impulse,  Is, 
or  exhaust  velocity,  v.  Specific  impulse  and  exhaust  velocity  are  related  via  the  equation, 
v  =  Is  go- 


30 


The  article  notes  that  several  different  organizations  have  conducted 
optimizations  for  various  missions  and  that  those  optimizations  assumed  a  constant 
specific  power  of  the  propulsion  system.  This  assumption  is  arguably  a  potential  source 
of  error  in  that  it  implies  the  following  about  variations  in  power  levels  and  specific 
impulse.  First,  it  assumes  that  constant  thrust  and  system  efficiency  could  be  evident 
with  a  changing  specific  impulse.  Additionally,  it  implies  that  constant  specific  power 
could  be  maintained  with  changing  power  levels  during  optimization.  (Auweter-Kurtz  et. 
al.,  1988,  pg.  512) 

In  undertaking  this  optimization  analysis,  this  point  was  certainly  considered. 
However,  due  to  the  fact  that  the  trace  of  an  engine’s  efficiency  is  not  continuous, 
average  values  were  utilized  over  limited  specific  impulse  ranges  thereby  assuming  no 
strong  dependence  on  specific  impulse  or  exhaust  velocity. 

Still,  it  is  not  elementary  to  assign  a  values  to  a  particular  engine. 

Items  of  concern  in  determining  specific  power  include  such  things  as  the  number 
of  engines  that  will  be  necessary  for  a  specific  mission  and  whether  or  not  those  same 
engines  can  share  various  components  of  the  power  conditioning  equipment.  If  power 
conditioning  components  can  be  shared,  then  a  improves  as  mass  is  decreased. 

Likewise,  a  would  be  detrimentally  affected  if  more  mass  were  required  such  as 
additional  solar  cell  panels  to  power  required  thrusters. 

Specific  power  also  can  be  improved  if  the  electrical  propulsion  system  can  be 
powered  through  the  systems  already  integral  to  the  payload.  As  such,  the  question  arises 
as  to  where  that  mass  should  be  billed.  If  it  is  attributed  to  the  payload  mass  (mpi),  then 
again  the  value  for  specific  power  is  improved. 
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In  summary  of  the  above  concerns,  it  should  be  evident  that  the  a  parameter  is  not 
only  engine  specific  but  also  configuration  (number  of  engines  and  shared  components  or 
redundancy)  and  design  (payload  compatible  power  sources)  specific. 

In  order  to  demonstrate  this  aspect  of  the  determination  of  the  value  for  specific 
power,  the  following  example  is  illustrated. 

Atlantic  Research  Corporation  is  one  of  the  current  suppliers  of  one  of  the 
commercially  available  Hall-type  thrusters.  In  their  brochure  for  the  SPT-100,  they 
depict  the  configuration  for  a  “typical  SPT-100  propulsion  system  for  GEO  satellites.” 

The  system  depicted  consists  of  the  following  major  components:  four  SPT-100 
thrusters,  four  pairs  of  xenon  flow  controllers  (XFC’s),  two  PPU-100  power  processing 
units,  the  propellant  management  assembly  (PMA),  and  two  xenon  storage  tanks.  The 
entire  configuration  allows  for  redundancy;  therefore,  we  would  assume  that  the  designed 
mission  would  require  the  use  of  only  two  of  the  four  thrusters.  (ARC  brochure,  1998) 

In  determining  specific  power,  it  is  necessary  to  account  for  the  masses  of  all 
components  that  encompass  the  power-plant  mass  (mpp).  This  does  not  include  the  actual 
mass  of  the  propellant  (mp).  Thus,  it  is  often  referred  to  as  the  “dry  mass”  --  probably 
hailing  from  the  early  days  of  liquid  rocket  propulsion. 

The  masses  of  all  of  the  components  are  as  follows: 

SPT-100  (4) . 3.5  kg  /  SPT 

XFC  (4  pairs) . 0.64  kg  /  pair 

PPU-100  (2) . 6.2  kg  /  PPU 

PMA . . 3.25  kg 

However,  the  mass  of  the  xenon  tanks  is  mission  dependent.  (Clauss,  1999) 
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The  tank  mass  can  be  estimated  in  the  manner  demonstrated  in  Chapter  II,  C. 
Again,  assuming  that  a  tank  which  holds  zero  propellant  mass  weighs  1  kg  and  that  the 
tank  has  a  mass  fraction  of  0.1  to  0.12  (as  a  function  of  the  mass  of  xenon),  then  for  a  40 
kg  of  xenon,  a  5  kg  tank  should  be  required  as  1  kg  +  (0.1)(40  kg)  =  5  kg.  (Clauss,  1999) 

In  the  brochure  example  for  a  GEO  satellite,  we  pick  a  usable  example  of  a  15- 
year  GEO  satellite  for  NSSK.  Such  a  mission  would  require  approximately  120  kg  of 
xenon  to  be  held  in  two  tanks.  Hence,  60  kg  of  xenon  would  be  held  in  each  tank. 
Accordingly,  one  tank  would  then  weigh  7  kg  as  1  kg  +  (0.1)(60  kg)  =  7  kg. 

Nowhere  in  any  current  publications  is  there  a  discussion  as  to  how  to  “charge” 
mass  values  in  determining  specific  power  or  “alpha.”  In  this  case,  all  “dead”  mass  is 
attributed  to  power-plant  mass  (mpp),  and  all  mass  that  is  proportional  to  power  is 
assigned  to  payload  mass  (mpi).  Hence,  we  would  have  to  include  2  kg  of  tankage  in  our 
value  for  power-plant  mass. 

With  the  value  for  tank  mass,  it  is  now  possible  to  calculate  the  power-plant  mass 
as  follows: 

mpp  =  4  (3.5  kg)spT  +  4  (0.64  kg)xFC  +  2  (6.2  kg)ppu  +  3.25  kg  pma  +  2.0  kg  tanks- 
Thus,  mpp  =  34.2  kg. 

Each  of  the  SPT-100  thrusters  requires  1 .35  kW  of  power  delivered  at  300  V  and 
4.5  A.  With  a  known  PPU  efficiency  of  r|Ppu  =  0.93,  we  are  able  to  calculate  that  a  total 
of  5,808  W  must  be  delivered  to  the  thruster  bank.  Hence,  we  are  now  able  to  calculate 
the  specific  power  for  the  configuration.  (Clauss,  1999) 

a  =  Pe  /  mpp  =  5,808  W  /  34.2  kg  =  169.8  W/kg. 
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It  can  not  be  overemphasized  that  this  calculated  value  of  specific  power  for  the 
SPT-100  is  only  for  the  illustrated  four-thruster  configuration  as  designed  for  the  mission 
of  a  15-year  GEO  satellite.  Otherwise,  the  masses  change,  and  specific  power  must  be 
re-calculated. 
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Table  5:  Summary  of  Current  Technology  Engines 
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IV.  PLANNING  /  OPTIMIZATION  ALGORITHM  -  MATCHING  THE 


MISSION,  PAYLOAD,  AND  ENGINE 


OPTIMUM  DESIGN  PROCEDURE: 


1.  Specify  the  Mission. 

Define  the  mission  profile  in  terms  of  the  velocity  requirements  (Au)  and  payload 
mass  (mpi).  List  additional  constraints  such  as  mission  time  (tm),  power  (Pe),  cost,  etc. 

2.  Estimate  the  payload  mass  ratio,  mpi  /  m0. 

Based  upon  past  operational  and  flight-test  experience,  an  initial  estimate  of 
suitable  candidate  electric  propulsion  thrusters  may  yield  an  approximate  expected 
payload  mass  ratio.  Note:  This  initial  “ballpark”  estimate  is  primarily  based  upon 
mission  profile  (velocity  requirement)  and  suitable  engines,  looking  at  cost,  mission  time, 
and  operational  flight  experience. 

In  order  to  remain  near  the  optimum  for  engine  performance  (as  illustrated  in 

earlier  plots),  a  potential  payload  mass  ratio  range  (from  mpi  =  0.35  to  0.55)  is  selected 

with  the  corresponding  numbers  for  Au  /  vc  and  v  /  vc. 

IDEAL:  without  “tankage”  penalty  (i.e.,  particularly  for  Teflon  PPT): 

mpi  /  m0  =  0.35:  Au  /  vc  =  0.404  ;  v  /  vc  =  0.782. 

mpi  /  rrio  =  0.45:  Au  /  vc  =  0.327  ;  v  /  vc  =  0.820. 

mpi  /mo=  0.55:  Au  /  vc  =  0.257  ;  v  /  vc  =  0.864. 

“TANKAGE”:  with  10  %  tankage  penalty: 

Au  /  vc  =  0.375  ;  v  /  vc  =  0.850. 

Au  /  vc  =  0.306  ;  v  /  vc  =  0.895. 

Au/vc  =  0.240;  v/vc  =  0.945. 


mpi  /  nio  =  0.35 
mpi  /  mo  =  0.45 
mpi  /  m0  =  0.55 
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If  none  of  the  listed  payload  ratios  will  suffice  (or  if  further  iteration  shows  that 
an  intermediate  value  for  payload  ratio  is  appropriate),  the  required  numbers  for  Au  /  vc 
and  v  /  vc  can  be  obtained  from  the  plots  in  Figures  8  and  9.  These  graphs  are  for  ideal 
systems  and  systems  requiring  “tankage”  respectively. 

3.  Calculate  the  corresponding  specific  impulse,  Is,  based  upon  payload  mass 
ratio. . .  using  the  tabular  data  below  and  the  equation, 

Is  =  V  /  go. 

IDEAL:  without  “tankage”  penalty: 

mpi /mo  =  0.35:  Is  =  1.936  Au/ g0. 

mpi  /  mo  =  0.45:  Is  =  2.5 1  Au  /  g0. 

mpi  /  mo  =  0.55 :  Is  =  3.36  Au  /  g0,  where  g0  =  9.8 1  m  /  s2. 

“TANKAGE”:  with  10  %  tankage  penalty: 

mp]  /  nio  =  0.35:  Is  =  2.27  Au  /  g0. 

mpi  /  mo  =  0.45 :  Is  =  2.92  Au  /  g0. 

mp|/mo  =  0.55 :  Is  =  3.94  Au / g0,  where  g0  =  9.8 1  m / s2. 

For  intermediate  values  of  payload  mass  ratio,  this  same  calculation  can  be 
accomplished  using  the  ratios  from  Figures  8  and  9  along  with  the  following  equation, 

Is  =  (  v  /  vc )  (  vc  /  Au  )  (  Au  /  go  ). 

4.  Designate  a  candidate  electric  propulsion  engine. 

Referencing  Figure  8  (on  page  35),  select  the  engine  (or  engines)  whose  available 
range  of  specific  impulse,  Is,  most  closely  matches  the  calculated  value  necessary  for  the 
optimum  profile. 

Engine  selection  also  locks  in  the  characteristic  engine  performance  values  for 
specific  power  (a),  thruster  efficiency  (r|t),  and  available  thrust.  These  values  can  be 
found  in  the  comprehensive  summary  of  Table  5  on  pages  36  and  37. 
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Again,  if  no  engine  proves  an  adequate  match  for  the  calculated  values  of  Is,  it  is 
necessary  to  use  the  graphs  depicted  in  Figures  9  and  10  and  attempt  intermediate  values 
of  mass  payload  ratio  within  the  range. 

5.  Calculate  burn  time,  tp. 

Determine  the  thruster  bum  time  for  the  profile  from  the  following  formula, 
tP  =  (( vc  /  Au )  Au)2  /(2  a  r|t). 

At  this  point  in  the  process,  all  of  the  above  values  are  given:  Au  from  mission 
profile,  vc  /  Au  from  step  2,  a  and  r|t  from  Table  5. 

6.  Check  that  burn  time,  tp,  is  less  than  both  mission  time,  tm,  as  well  as  the  design 
lifetime  of  the  designated  engine. 

If  bum  time  exceeds  acceptable  mission  time  or  design  life  of  the  thruster,  select 
another  electric  propulsion  candidate  (step  4)  or  vary  the  payload  mass  ratio  (step  2). 

7.  Calculate  the  anticipated  masses  of  the  propellant  (mp)  and  power-plant  (mpp)  via 
the  following  formulas  for  the  respective  payload  mass  ratios: 

IDEAL:  without  “tankage”  penalty  (i.e.,  for  Teflon  PPT): 

mpi  /  mo  =  0.35:  mp  =  1.152  mpi ;  mpp  =  0.612  mp 

mpi  /  nio  =  0.45:  mp  =  0.73 1  mpi ;  mpp  =  0.672  mp 

mpi  /  mo  =  0.55:  mp  =  0.469  mp] mpp  =  0.746  mp 

“TANKAGE”:  with  10  %  tankage  penalty  (i.e.,  for  all  other  engines): 

mpi  /  mo  =  0.35:  mp  =  1.078  mpi ;  mpp  =  0.723  mp 

mpi  /  nio  =  0.45:  mp  =  0.679  mpi ;  mpp  =  0.801  mp 

mpi  /  nio  =  0.55:  mp  =  0.432  mP] ;  mpp  =  0.893  mp 

For  intermediate  payload  ratios,  the  following  equations  must  be  used... 

(a)  mpi  /  mo  =  desired  ratio. 

(b)  mpp  =  ( v  /  vc  )2  mp. 

(c)  mo  =  mp+  mpp  +  mpi. 

(d)  mpi  is  given. 
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8.  Determine  the  total  mass  of  the  assembled  vehicle  (m0)  and  verify  the  payload 
mass  ratio. 

mo  =  mp  +  mpp  +  mpi.  mass  ratio  =  mpi  /  m0. 

9.  Determine  the  required  thrust  based  on  the  values  for  specific  impulse,  Is,  and 
propellant  mass,  mp.  Determine  the  required  power  to  be  delivered  to  the  thruster. 

F  =  ( mp  /  tp )  Is  go.  Pe  =  a  mpp. 

Note:  This  calculated  value  of  Pe  should  not  greatly  exceed  the  current  “state  of 
the  art”  for  power  generation  (fuel  cells,  battery,  or  solar  cells)  of  approximately  20  kW. 
If  the  power  requirements  cannot  be  met,  return  to  step  4. 

10.  From  the  thrust  calculation  and  the  engine  thrust  characteristics  (step  4), 
calculate  the  number  of  engines  required  to  provide  that  thrust  for  the  assigned 
mission. 

11.  Design  check. 

Considering  component  redundancy,  is  this  design  acceptable?  That  is,  with 
current  technology,  is  this  a  reasonable  number  of  required  engines?  If  not,  return  to 
design  step  4  or  attempt  an  intermediate  payload  ratio  (step  2). 

If  yes,  it  may  be  advisable  to  iterate  in  order  to  account  for  redundancy  thereby 
improving  the  a  parameter.  In  other  words,  are  there  any  components  that  the  thrusters 
may  be  able  to  share,  thereby  reducing  mpp  and  improving  a.  This  will  improve  overall 
performance  by  reducing  mass  and  possibly  reducing  the  number  of  engines  required  for 
the  mission. 
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0  0.3  0.34  0.36  0.38  0.4  0.42  0.44  0.46  0.43  OvS  0.52  0.54  O.SS  G.53  0.6 

tripi  /  m0 

Figure  9:  Comprehensive  Plot  of  Au  /  vc  and  v  /  vc  for  Ideal  Systems. 

This  plot  serves  as  a  baseline  plot  of  the  values  of  Au  /  vc  and  v  /  vc  for  all 
anticipated  payload  mass  ratios  between  the  optimum  selection  points  of 
mpi  /  mo  =  0.35  and  0.55.  The  points  roughly  bracket  the  optimum  point  at 
mpi  /  m0  =  0.45. 

This  plot  is  usable  only  for  ideal  estimates  of  the  calculations  discussed  in  the 
design  calculation,  and  it  is  usable  for  systems  that  do  not  require  fuel  tanks  -  i.e.  no 
“tankage  penalty.”  Such  electric  thrusters  or  engines  include  die  pulsed  MPD  thrusters: 
Teflon  PPT,  LES  8/9  PPT,  NASA  Primex  EO-1,  and  the  Primex  PRS-101. 
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Figure  10:  Comprehensive  Plot  of  Au  /  vc  and  v  /  vc 
for  Systems  Requiring  “Tankage”  of  10%. 


This  plot  is  acceptable  for  use  with  all  electric  propulsion  thrusters  that  require 
fuel  storage  tanks.  Therefore,  it  may  be  used  for  “tankage  penalty”  calculations  for  all  of 
the  engines  listed  except  the  pulsed  MPD  type. 

Because  of  the  fact  that  these  plots  take  fuel' tank  considerations  into  account,  this 
chart  (or  at  least  the  determined  values  for  the  defined  payload  rations)  should  be 
consulted. 


44 


46 


V. 


APPLICATION  OF  METHOD 


Note:  In  the  following  three  examples,  the  defined  parameters  which  had  to  be  estimated 
by  the  author  are  preceded  by  the  symbol 
A.  “LIGHT”  LEO  TO  GEO  MISSION 

For  a  particular  LEO  to  GEO  satellite,  the  following  parameters  are  defined: 
mission  velocity  requirement,  Au  =  4200  m/s 
payload  mass,  mpi  =  250  kg 
mission  time,  tm  =  unrestricted  /  reasonable. 

Using  the  defined  optimization  procedure,  three  potentially  optimized  engine  runs 
are  attempted  at  payload  mass  ratios  of  mpi  /  mo  =  0.35, 0.45,  and  0.55. 

In  this  example,  the  second  (“tankage”  considered)  case  for  the  payload  mass 
ratio  of  0.45,  proves  to  be  the  optimized  profile.  Note  that  approximately  10  to  12 
engines  are  required  if  the  H2  arcjet  is  selected  as  the  thruster  for  this  profile. 

This  example  demonstrates  that  fuel  tank  mass  considerations  lead  to  the  addition 
of  more  thrusters  in  order  to  compensate  for  the  added  mass.  The  mass  penalty  is  easily 
scalable. 
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IDEAL: 


mpi  /  mo  =  0.35 

Is  =  829  sec 
engine  selected: 

NH3  Arcjet 

Is  =  650  sec 
a  =  320  W/kg 
rit  =  0.32 
T  =  0.2  N 

tp=5.28  x  105  sec 
=  6.11  days 
tp  <'  tm 

mp  =  270  kg 
mpp=  195.2  kg 
mo  =  715  kg 
payload  ratio  checks 

F  =  3.26  N;  Pe  =  63  kW 
17  engines  required; 
power  req.  excessive 

with  TANKAGE... 

Is  =  972  sec 
engine  selected: 

H2  Arcjet 

Is  =  1000  sec 
a  =  333 
T|t  =  0.4 

T  =  0.2  -  0.25  N 

tp=  4.71  x  105sec 
=  5.45  days 
tp  <-  tm 

nip  =  270  kg 
mpp  =195.2  kg 
mo  =  715  kg 
payload  ratio  checks 

F  =  5.47  N;  Pe  =  65  kW 
.*.  22  to  28  engines 
required 


nip]  /  m0=  0.45 

Is=  1075  sec 
engine  selected: 

NASA  /  Primex  EO-1 
Is=  1150  sec 
a  =  20  W/kg 
r(t  =  0.098 
T  =  pulsed 

tp  =  4.21  x  107  sec 
=  487  days 

tp  ^  tm 

mp=  182.8  kg 
mpp=  122.8  kg 
mo  =  556  kg 
payload  ratio  checks 

F  =  0.0490  N;  Pe  =  3kW 
.’.  1,919  engines 
required!!! 

Is  =  1250  sec 
engine  selected: 

H2  Arcjet 

Is  =  1000  sec 
a  =  333 
t|t  =  0.4 

T  =  0.2  -  0.25  N 

tp=  7.07  x  105  sec 
=  8.18  days 
tp  ^  tm 

mp  =  169.8  kg 
mpp=  136.0  kg 
mo  =  556  kg 
payload  ratio  checks 

F  =  2.36N;  Pe  =  46kW 
.*.  10  to  12  engines 
required;  power  high 


mpi  /  nio=  0.55 

Is  =1439  sec 
engine  selected: 

ARC  SPT-100 

Is  =  1600  sec 
a  =  ~  169.8W/kg 
Pt  =  0.48 
T  =  0.083  N 

tp=  1.638  x  106  sec 
=  19.0  days 
tp<tm 

mp=  117.3  kg 
mpp  =  87.5  kg 
mo  =  455  kg 
payload  ratio  checks 

F=  1.011  N;  Pe=15kW 
.•.  13  engines  required 


Is=  1687  sec 
engine  selected: 

Hall  (Xe) 

Is  =  15-1600  sec 
a  =  150 
rit  =  0.5 
T  =  0.04  N 

tp=2.04  x  106  sec 
=  23.6  days 

tp  <-  tm 

mp=  108.0  kg 
mpp  =  96.4  kg 
mo  =  454  kg 
payload  ratio  checks 

F  =  0.83  IN;  Pe=  15  kW 

.*.  21  engines  required 
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B.  “HEAVY”  LEO  TO  GEO  (COMM  SAT)  MISSION 

Given  a  LEO  to  GEO  mission  for  the  orbit  raising  of  a  communication  satellite, 
the  optimization  algorithm  requires  that  we  initially  define  the  payload  and  required  Au 
for  the  transit. 

The  payload  is  defined  to  be  relatively  massive,  weighing  1000  kg,  and  the 
defined  Au  for  the  transfer  is  approximately  4640  m/s.  For  this  mission,  the  transfer  time 
cannot  be  unrealistically  excessive  due  to  the  general  mission  urgency  of  communication 
satellite  tasking.  Using  the  planning  algorithm,  the  results  are  illustrated  as  shown  on  the 
following  page. 

This  example  provides  some  interesting  results.  First,  the  payload  is  extremely 
heavy  for  most  electrical  propulsion  engines.  As  stated  in  the  earlier  chapters,  most  high 
thrust  electric  propulsion  engines  are  still  in  development.  Yet,  the  Hall  TAL  D-55 
engine  is  capable  of  accomplishing  this  mission  in  88)4  days  with  only  1 1  engines  - 
excluding  redundancy. 

In  both  of  the  cases  at  the  payload  mass  ratio  of  0.45,  it  is  evident  that  the 
increased  bum  time  to  maintain  the  profile  allows  for  fewer  engines  as  less  thrust  is 
required  over  a  longer  time.  None-the-less,  the  number  of  engines  required  for  this 
profile  as  mpi  /  mo  =  0.45  is  unacceptable. 
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IDEAL: 


mpi  /  mo  =  0.35 

Is  =  916  sec 

engine  selected: 

NASA/Primex  EO-1 
Is  =  1150  sec 
a  =  ~  20  W/kg 
ri,  =  0.098 
T  =  pulsed 

tp=  3.37  x  107sec 
=  389  days 
=  1.067  years 

SHOWSTOPPER... 
due  tp » tm 


with  TANKAGE... 

Is  =  1074  sec 
engine  selected: 

H2  Arcjet 

Is  =  1000  sec 
a  =  333  W/kg 
rit  =  0.4 

T  =  0.2  -  0.25  N 

tp=  5.75  x  105  sec 
=  6.65  days 

tp  tm 

nip  =  1078  kg 
mpp  =  779  kg 
mo  =  2857  kg 
payload  ratio  checks 

F=  18.39  N;Pe=  259  kW 
/.  74  to  92  engines 
required! ! ! 


mpi  /  mo  =  0.45 

Is=  1187  sec 
engine  selected: 

Hall  (Xe) 

Is  =15 -1600  sec 
a  =150  W/kg 
rit  =  0.5 
T  =  0.04  N 

tp=  1.342  x  106  sec 
=  15.54  days 

tp  <'  tm 

mp  =  731  kg 
mpp  =  491  kg 
m0  =  2222  kg 
payload  ratio  checks 

F  =  8.02  N;  Pe  =  74kW 
/.201  engines  required! 


Is=  1381  sec 
engine  selected: 

Hall  (Xe) 

Is=  15  -1600  sec 
a  =  150  W/kg 
rit  =  0.5 
T  =  0.04  N 

tp=  1.533  xlO6  sec 
=  17.74  days 

tp  <'  tm 

mp  =  679  kg 
mpp  =  544  kg 
mo  =  2223  kg 
payload  ratio  checks 

F  =  6.52  N;  Pe  =  82  kW 
.’.  163  engines  required! 


mpi  /  nio=  0.55 

Is  =  1589  sec 
engine  selected: 

Hall  TAL  D-55 

Is  =  1600  sec 
a  =  ~  50.9  W/kg 
Tit  =  0.48 
T  =  0.082  N 

tp=  6.67  x  106  sec 
=  77.1  days 

tp  ^  tm 

mp  =  469  kg 
mpp  =  350  kg 
mo=  1819kg 
payload  ratio  checks 

F  =  1.104  N;  Pe=18kW 
14  engines  required 


Is=  1864  sec 
engine  selected: 

Hall  TAL  D-55 

Is=  1600  sec 
a  =  ~  50.9  W/kg 
Tit  =0.48 
T  =  0.082  N 

tp=  7.65  x  106  sec 
=  88.5  days 

tp  ^  tm 

mp  =  432  kg 
mpp  386  kg 
mo  =  1818  kg 
payload  ratio  checks 

F  =  0.886  N;  Pe=  20  kW 
.*.  1 1  engines  required 
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c. 


LEO  TO  MARS  MISSION 


For  a  representative  LEO  to  Mars  mission,  the  optimization  algorithm  requires 
that  we  initially  define  the  payload  and  required  Au. 

The  payload  is  defined  to  be  a  relatively  small  probe  (mass  equal  to  250  kg),  and 
the  defined  Au  for  the  transfer  is  approximately  5,700  m/s.  In  this  case,  Au  also  defines  a 
mission  time  (tm)  equal  to  256  days.  Thus,  in  order  for  the  calculations  to  be  valid,  the 
bum  time  (tb)  must  not  exceed  tm. 

The  results  of  the  calculations  are  demonstrated  as  shown. 

Note  that  for  the  third  case  of  mpi  /  mo  =  0.55  either  of  the  steady  MPD  thrusters 
would  have  been  a  better  specific  impulse  fit.  Those  engines  are  scaleable  in  the  range 
from  Is  =  2000  to  5000  sec.  However,  since  they  are  presently  still  under  development, 
they  have  been  excluded  from  this  analysis  even  though  it  is  strongly  anticipated  that 
such  engines  will  be  the  desired  propulsion  systems  for  deep-space  missions. 
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IDEAL: 


mpi  /  rrio=  0.35 

Is  =1125  sec 

engine  selected: 

Teflon  PPT  (pulsed) 
Is  =  1000  sec 
a  =  1.0  W/.kg 
lit  =  0.07 
T  =  pulsed 

tp=  1.422  x  109sec 
=  16,457  days 
=  45  years!!! 

SHOWSTOPPER... 
due  tp  » tm 


with  TANKAGE... 

Is  =  1319  sec 
engine  selected: 

Hall  (Xe) 

Is  =  15-1600  sec 
<x  =  150  W/kg 
Tit  =  0.5 
T  =  0.04  N 

tp=  1.540  x  106  sec 
=  17.83  days 

tp  <  tm 

mp  =  270  kg 
mpp=  195.2kg 
mo  =  715  kg 
payload  ratio  checks 

F  =  2.58  N;  Pe  =  30kW 
/.  65  engines  required; 
power  required  is  high 


mpi  /  mo  =  0.45 

Is  =  1458  sec 
engine  selected: 

Hall  SPT  (Xe) 

Is=  1500  sec 
a  =150  W/kg 
r]t  =  0.5 
T  =  0.04  N 

tp  =  2.03  x  106  sec 
=  23.4  days 

tp  <  tm 

mp  =  182.7  kg 
mpp  =  122.8  kg 
mo  =  556  kg 
payload  ratio  checks 

F=  1.324  N;  Pe=19kW 
.’.  34  engines  required 


Is  =  1697  sec 
engine  selected: 

ARC  SPT- 100 

Is  =  1600  sec 
a  =  131.4  W/kg 
lit  =  0.48 
T  =  0.083  N 

tp  =  2.75  x  106  sec 
=  31.8  days 
tp  tm 

mp  =  169.8  kg 
mpp=  136.0  kg 
nio  =  556  kg 
payload  ratio  checks 

F  =  0.969  N;Pe=  18  kW 
.’.  12  engines  required 


mpi  /  m0=  0.55 

Is=  1952  sec 
engine  selected: 

SPT- 100  (Fakel) 

Is  =  1600  sec 
a  =150  W/kg 
lit  =  0.48 
T  =  0.083  N 

tp=  3.42  x  106  sec 
=  39.5  days 

tp  ^  tm 

mp=  117.3  kg 
mpp  =  87.5  kg 
mo  =  455  kg 
payload  ratio  checks 

F  =  0.538  N;Pe=  14  kW 
.*.  7  engines  required 


Is  =  2289  sec 
engine  selected: 

XIPS 

Is  =  28-3500  sec 
a  =  100  W/kg 
Tit  =  0.75 
T  =  0.015-0.04  N 

tp=  3.76  x  106  sec 
=  43.5  days 

tp <-  tm 

mp  =  108.0  kg 
mpp  =  96.4  kg 
mo  =  454  kg 
payload  ratio  checks 

F  =  0.789  N;  Pe=10kW 
/.20  to  53  engines 
required 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Routine  implementation  of  electric  propulsion  for  space  thrusting  applications  is 
relatively  new  even  though  the  governing  concepts  are  almost  40  years  old.  The  very 
latest  engines,  presented  as  examples  in  this  thesis,  demonstrate  that  some  aspects  of 
electric  propulsion  are  still  in  their  early  stages  of  development.  However,  in  light  of 
time-tested  theory  and  the  limitations  of  chemical  propulsion  systems,  it  is  evident  that, 
right  now,  electric  propulsion  is  the  only  potential  stepping-stone  to  deep-space  probes  or 
other  craft  that  could  conceivably  travel  the  great  distances  between  the  planets  and 
beyond  the  solar  system. 

For  any  such  mission  profiles  -  interplanetary  or  interstellar  in  particular  -  the 
propulsion  systems  for  these  missions  will  have  to  be  optimized.  Such  engines  will  have 
to  be  cost-effective,  thrust  efficient  and  therefore  mass  minimized.  Someday  there  may 
even  be  considerations  to  make  these  systems  potentially  re-useable  and  refuelable.  This 
thesis  presents  one  technique  with  which  mission  profiles  and  engine  selection  can  be 
matched  and  consequently  optimized  to  theoretically  achieve  the  greatest  benefit  of  least 
overall  mass  with  both  maximum  fuel  and  maximum  payload  in  order  to  achieve  the 
shortest  mission  time  for  a  given  payload. 

Thus,  a  primary  focus  of  this  thesis  has  been  to  provide  a  new  way  of  looking  at 
electrical  propulsion  systems  design  in  addition  to  attempting  to  provide  an  up-to-date 
and  comprehensive  summary  of  available  electric  propulsion  thruster  data.  The  selected 
approach  (i.e.,  the  Langmuir-Irving  formulation)  is  based  upon  first  principles  as  well  as 
an  analysis  of  other  methods  that  are  currently  in  use.  In  essence,  an  objective  has  been 
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to  produce  a  design  alternative  to  the  satellite  design  engineers.  This  methodology  or 
“procedure”  is  intended  to  generate  some  consideration  into  performance  near  a  pre¬ 
determined  yet  not  inflexible  “dual-optimum”  point.  It  is  apparent  in  the  literature  that 
aspects  of  the  discussed  “dual-optimum”  point  have  been  plotted  in  various  texts  and 
somewhat  addressed  but  never  fully  explored. 

The  most  important  concept  to  note  is  the  idea  that  these  results  are  “theoretical” 
in  that  systems  are  currently  designed  according  to  algorithms  that  are  based  more  solidly 
on  operational  experience  and  power  considerations  than  on  any  approach  that  considers 
pure  optimization.  In  essence,  no  design  algorithm  is  really  in  use  that  seems  to  consider 
optimization.  Granted,  electric  propulsion  is  still  in  its  early  stages  and  most  applications 
involving  satellite  re-positioning  and  orbit-lifting  do  not  place  heavy  lift  nor  aggressive 
profile  constraints  upon  engines.  Consequently,  such  considerations  or  additional 
restraints  may  not  really  be  required. 

Due  to  the  relatively  limited  number  of  electrical  propulsion  engines  that  are 
currently  operational,  the  examples  of  this  thesis  show  that  it  is  sometimes  difficult  to 
find  an  identical  match  to  an  optimum  profile’s  specific  impulse,  Is.  In  fact,  some  of  the 
most  powerful  and  most  promising  electric  propulsion  engines  -  those  most  likely  to  be 
desirable  for  high  thrust  and  high  payload  mass  missions  -  are  still  in  some  stages  of 
development. 

Specific  impulse,  Is,  is  comparable  to  the  required  “gas  mileage”  for  a  particular 
trip  through  space.  After  the  profile  is  determined  through  orbital  mechanics,  the 
required  incremental  change  in  vehicular  velocity  (Au)  is  known.  Then,  the  required  tank 
of  fuel  can  be  established  for  an  expected  level  of  performance  or  fuel  efficiency  -  “gas 
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mileage.”  Hence,  the  objective  is  to  match  the  profile  to  specific  impulse  (or  engine). 
This  is  accomplished  using  the  “dual-optimum”  procedure.  Nonetheless,  the  “dual¬ 
optimum”  establishes  not  only  Is  but  also  suggests  desired  values  of  specific  power  (a) 
and  thruster  efficiency  (pt)  as  well  as  the  likely  bum  time  (tp)  for  a  desired  mission 
profile,  Au. 

As  demonstrated  in  the  three  provided  examples  in  the  previous  chapter,  it  is  not 
obvious  how  a  particular  engine  will  perform  on  a  specific  mission  profile.  In  essence, 
the  design  optimum  is  dependent  as  much  on  the  payload  mass  (mpi)  and  mass  payload 
ratio  (mpi  /  nio)  as  upon  the  specific  characteristic  of  the  engine  itself. 

The  most  important  of  these  engine  characteristics  are  specific  power  (a),  thruster 
efficiency  (r|t),  and  available  thrust  (T).  Specific  power,  a,  is  the  most  critical  to  the 
designer  because  it  drives  the  amount  of  available  power  to  the  thruster  versus  the 
required  mass  of  the  overall  system  that  provides  that  power.  As  power  conditioning 
systems  improve  and  the  “state-of-the-art”  in  the  field  advances,  this  parameter  will 
improve  and  enhance  electrical  propulsion  system  performance.  In  this  thesis,  it  is 
evident  that  specific  power  has  been  the  most  difficult  parameter  to  establish  for  the 
majority  of  the  listed  engines.  A  worked  example  demonstrates  that  this  is  primarily  due 
to  the  fact  that  the  value  given  to  a  varies  depending  somewhat  upon  arbitrary 
considerations  such  as  thruster  configuration,  redundancy,  etc.  In  addition,  it  shows  that 
mass  reduction  is  equally  as  important  as  where  that  mass  is  “billed”  within  the 
calculations. 

Thruster  efficiency,  r|t ,  generates  a  feel  for  the  engine’s  efficiency  in  terms  of 
nozzle  losses,  off-axis  losses,  as  well  as  other  thruster  inefficiencies.  Generally,  these 
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values  vary  widely  depending  upon  the  type  of  thruster  under  consideration.  For 
example,  the  LES  8/9  PPT  thruster  has  an  efficiency  of  only  0.0068  to  0.009  while  the 
N2H4  resistojet  turns  in  numbers  in  the  0.8  to  0.9  range.  The  reasons  for  this  are 
discussed  in  depth. 

Available  thrust,  T,  although  extremely  small  on  most  electrical  propulsion 
systems,  gives  a  measure  of  the  number  of  engines  required  for  a  necessary  profile.  This 
value  can  be  a  major  player  as  demonstrated  in  the  examples  because  it  drives  largely  the 
additional  issues  of  component  redundancy  and  additional  mass  -  which  again  affects 
alpha,  a. 

Another  consideration  which  became  apparent  in  the  course  of  this  work  is  the 
capability  to  allow  for  the  required  fuel  tanks  within  the  design  and  thus  within  the 
optimum  profile.  Fuel  considerations  are  certainly  mission  or  profile  dependent. 
Consequently,  this  is  a  substantial  capability  for  the  designer  to  predict  optimum 
performance  with  variable  tank  mass  as  per  mission  requirements. 

These  results,  as  presented,  are  easily  scalable  depending  upon  the  particular 
profile  and  its  fuel  requirements.  Exploration  of  the  topic  of  allowing  for  “tankage” 
created  estimation  possibilities  that  allowed  for  even  more  realistic  comparisons  between 
different  electric  propulsion  systems.  As  evident  in  the  overall  engine  parameter 
summary  in  Table  5,  some  engines  (of  the  pulsed  MPD  variety)  require  no  fuel  storage 
tanks  at  all. 

If  the  engine  summary  shows  one  thing,  it  certainly  lets  one  know  that  there  are  a 
host  of  possibilities  out  there  being  developed  by  aggressive  companies.  These  systems 
use  a  variety  of  physical  principles  to  create  thrust.  Some  at  greater  expense  of  the  other 
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relevant  parameters  —  a,  r)t,  and  T.  Consequently,  the  examples  provided  for  the  design 
algorithm  (on  the  three  different  mission  profiles)  demonstrate  that  all  engines  are  not 
created  equal  for  various  mission  profiles  -  certainly  not  if  one  endeavors  to  fly  an 
optimum  profile. 

Of  course,  industry  may  argue  that  there  are  important  additional  considerations 
on  any  profile.  For  example,  a  communication  satellite  that  takes  six  months  to  be  in 
position  is  probably  not  cost-effective.  That  could  be  a  mission  failure.  Nonetheless, 
success  depends  on  what  parameter  is  being  optimized.  In  the  case  of  this  research,  the 
objective  was  to  optimize  both  the  mass  and  time  constraints.  At  the  present  time,  in  the 
development  of  the  “state-of-the-art”  in  the  evolution  of  electrical  propulsion,  this  was 
not  always  possible.  This  fact  is  certainly  demonstrated  in  the  examples  provided,  which 
show  that  there  are  no  obvious  trends  with  the  systems  currently  available.  Electric 
propulsion  is  not  totally  there  yet;  however,  scientists  and  engineers  are  pushing  the 
frontier  forward  and  expanding  to  new  horizons  and  possibilities. 

B.  RECOMMENDATIONS 

In  any  branch  of  research  involving  a  rapidly  growing  technology  (particularly 
one  in  which  corporate  vendors  on  the  “cutting-edge”  of  that  technology  are  especially 
unwilling  to  give  up  data  that  they  consider  proprietary),  it  is  extremely  difficult  to  track 
the  leading  edge  of  developments.  Consequently,  acquiring  the  latest  available  data  can 
be  next  to  impossible.  Hence,  best-estimate  scenarios  can  be  flawed  or  only  approximate. 

In  the  course  of  this  research,  the  author  frequently  ran  into  barriers  when 
attempting  to  answer  questions  regarding  particular  thruster  systems.  A  more  open  forum 
could  certainly  benefit  not  only  science  but  also  the  vendors  themselves,  particularly  if 
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the  data  contained  within  this  thesis  proves  useful  to  designers.  On  the  other  hand,  the 
developed  procedure  might  demonstrate  that  certain  engines  are  less  effective  for 
necessary  mission  profiles.  As  a  result,  one  design  might  prove  less  effective  than 
another. 

Additional  work  could  be  extremely  advantageous  if  a  more  in-depth  study  could 
be  made  into  attempting  to  combine  the  various  schools  of  thought  in  design  criteria. 

Like  missile  design,  there  is  no  clear-cut  or  established  method  that  is  followed  due  to  the 
fact  that  a  design  must  consider  the  whole  picture  -  i.e.,  mass,  time,  and  power 
considerations  -  all  at  once.  Considerations  of  power,  time,  and  cost  seem  to  be  the 
biggest  factors  in  design  at  the  present  time,  and  obviously  contractors  wanting  to  sell 
their  products  are  not  going  to  be  quick  to  admit  that  a  system  could  be  better  optimized 
with  another  company’s  thruster.  Certainly,  perfection  is  infinitely  expensive,  but  that  is 
not  what  this  process  advocates.  Systems  could  simply  be  better  if  any  attempt  were 
made  to  “ride  the  dual-optimum  profile.”  The  author  would  certainly  be  curious  if 
industry  were  to  consider  the  presented  ideas  and  improve  upon  them. 
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APPENDIX  A.  DEVELOPMENT  -  PAYLOAD  MASS  FRACTION  EQUATION: 


By  definition,  specific  power,  a  =  Pe  /  mpp,  Thus,  mpp  =  Pe  /  a.  (Sutton,  1992,  pg.  597) 

By  definition,  electrical  power  available  (as  previously  shown),  Pe  =  mp  v2  /  (2  tp  tp). 
Hence,  mpp  =  (1/a)  (mp  v2  /  (2  tp  r|t)).  Solving  for  payload  mass,  mp,  yields... 
mp  =  mpp  (2  a  tp  rit  /  v2)  or  mpp  /  mp  =  v2  /  2  a  tp  rjt.  (Sutton,  1992,  pg.  596) 

By  definition  of  allocated  spacecraft  masses,  mo  =  mp  +  mpi  +  mpp.  Define  final  mass,  mf 
=  mo  -  mp  (i.e.  initial  mass  less  propellant  mass).  Thus,  redefine  initial  mass  as  follows: 
mo  =  mp  +  mf.  (Sutton,  1992,  pg.  596) 

Derived  from  the  equation  for  maximum  velocity  at  propellant  burnout, 

Au  =  v  In  (mo  /  mf),  know  that  eAu/v  =  nio/  mf = m0  /  (m0  -  mp)  =  1  /  (1  -  mp/  mo).  This 
equation,  when  solved  for  mp  /  mo,  gives  mp  /  mo  =  1  -  e'Au/v.  (Sutton,  1992,  pg.  123) 

Rewriting  the  equation  for  initial  mass  yields  nio  =  mp  (1  +  mpp  /  mp)  +  mp|.  Now,  when 
substituting  the  above  equation  for  mpp  /  mp,  the  initial  mass  is  rewritten  as  follows: 

mo  =  mp  (1  +  v2  /  2  a  tp  rit )  +  mpi. 

Again,  rewriting  the  equation...  1  =  (mp  /  nio )  (1  +  v2  /  2  a  tp  t)t )  +  mpi  /  m0.  Substituting 
the  equation  for  mp  /  mo  and  solving  for  mpi  /  nio,  the  equation  transforms  to 
mp)  /  mo  =  1  -  (1  -  e"Au/v)  (1  +  v2  /  2  a  tp  r|t). 

When  this  equation  is  multiplied  through  and  solved  for  mpi/mo ,  the  payload  mass 
fraction  appears  as  illustrated:  mpi  /  nio  =  (l-(v/  vc)2(eAu/v  - 1))/  eAu/v. 
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APPENDIX  B.  FIRST  PRINCIPLES  -  CODES 


The  following  three  pages  contain  the  MATLAB  codes  that  were  utilized  in  the 
generation  of  Figures  2,  3,  and  4  respectively. 
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%  THESIS 

%  Plot2{2):  mpl  /mo  vs.  v  /  vc 

%  "payload  fraction  vs .  v  /  vc" 

%  LT  John  Jay  De  Beilis,  USN 
%  30  JUN  99  . 

%  variables: 

%  x ( )  =  v  /  vc 

%  a ( )  =  deltau  /  vc 

%  y()  =  mpl  /  mo 


xl=l inspace (0, 3, 300) 
x2=linspace (0,3,300) 
x3=linspace (0,3,300) 
al=0 . 2 
a2=0 . 3 
a3=0 . 5 

yl= ( 1-xl .  A2  . *  (exp (al . /xl ) -1) ) . /exp (al\  /xl ) 
y2= ( l-x2 . A2  . *  (exp (a2 . /x2) -I) ) . /exp (a2 . /x2 ) 
y3=(l-x3. A2  .*  (exp(a3./x3)-l) ) . /exp (a3 . /x3) 

figure (1) ,  plot (xl, yl, x2, y2,x3, y3) ,  grid 
axis ( [0  3  0.0  1.0]) 
ylabelt’mpl  /  mo 1 ) , xlabel ( 1 v  /  vc') 
gtext ( *deltau/vc  =  0.2* ),  gtext ( 1 deltau/vc  =0.3 


T),  gtext ( 1 deltau/vc  =  0.5') 
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%  THESIS 

%  Plotl:  deltau  /  vc  vs.  v  /  vc  (or  Is  /  vc) 

%  "incremental  change  of  vehicle  velocity  vs.  impulse" 

%  LT  John  Jay  De  Beilis,  USN 
%  25  AUG  99 

%  Variables : 

%  x()  =  v  /  vc 

%  a ( )  =  mpl  /  mo 

%  y()  =  deltau  /  vc 

xl=linspace (0, 3, 300) 
x2=linspace (0, 3, 300) 
x3=linspace (0, 3, 300) 

al=0 . 35 
a2-0 . 45 
a3=0 . 55 

yl=xl.*log(  (1  +  xl . A2 )  . /  (al  +  xl.A2)) 
y2=x2 . *log ( ( 1  +  x2.A2)./(a2  +  x2.A2)) 
y3=x3. *log ( (1  +  x3.A2)./(a3  +  x3.A2)) 

figure (1 ) ,  plot (xl, yl, x2, y2, x3, y3) ,  grid 

axis ([0.0  2.0  0.0  0.5] ) 

ylabel (' deltau  /  vc 1 ) , xlabel ( ' v  /  vc’) 

gtext ( ’mpl/mo  =  0.35f),  gtext ( 'mpl/mo  =  0.45’),  gtext ( ’mpl/mo  =  0.55* ) 
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THESIS 

Plot3:  optimization  curves  —  joint  optimum 

"product  of  incremental  change  of  vehicle  velocity  and  payload  ratio  vs, 

specific  impulse"  (mpl/mo  *  deltau  vs.  Is) 

LT  John  Jay  De  Beilis,  USN 
25  AUG  99 


Variables : 

x() 
a  ( ) 
y() 


=  V  /  VC 

=  mpl  /  mo 
=  deltau  /  vc 


xl=linspace (0,  3,  300) 
x2=linspace (0, 3, 300) 
x3=l inspace (0,3,300) 

al=0 . 35 
a2=0 . 45 
a3=0 . 55 

yl=xl . *log ( (1  +  xl . A2 ) .  / (al  +  xl.A2)) 
y2=x2 . *log ( ( 1  +  x2.A2)./(a2  +  x2.A2)) 
y3=x3 . *log ( ( 1  +  x3.A2)./(a3  +  x3.A2)) 

%zl=yl.*al 
%z2=y2 . *a2 
%z3=y3 . *a3 


figure (1) ,  plot (xl, yl, x2, y2, x3, y3) ,  grid 
axis ( [0.1  0.15  0.5  1.5] ) 

ylabel ( * deltau/vc  *  mpl/mo1),  xlabel ( 1 v/vc 1 ) 
gtext ( 1 mpl /mo  =  0.35')/  gtext ( ’mpl/mo  =  0.45’) 


gtext ( ’mpl/mo  =  0.55’) 
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APPENDIX  C.  DUAL  OPTIMUM  CODE 

This  appendix  lists  the  necessary  Maple  program  code  that  is  utilized  in  the 
production  of  Figure  4,  which  illustrates  the  “dual  optimum.” 
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THESIS:  Plot  #5:  OPT.  CURVES:  3-D  plot:  optimization  of  v/vc  and  mpl/mo  curves 

LT  John  Jay  De  Beilis,  USN 

30JUN  1999 

defined  variables: 

deltau  =  defined  incremental  change  of  vehicle  velocity  (m/s)  —  defined  by  mission 
vc  =  characteristic  velocity  (m/s) 
v  =  effective  exhaust  velocity  (m/s) 
mratio  =  payload  ratio  =  mpl/mo 

x  =  v/vc 
y  =  mpl/mo 

>  restart;  with(plots);  with(plottools) ; 

plot 3d  <x*y*ln ( ( 1+x*  2 ) / (y+x*  2)),x=0.5..1.3,y=0.3..0.6, axes=BOXED,  t i 
tle^optimization:  product  of  v/vc  and  mpl/mo"); 

[  animate,  animate3d,  animatecurve,  changecoords,  complexplot,  complexplot3d,  conformal, 

contourplot,  contourplot3d,  coordplot,  coordplot3d,  cylinderplot,  densityplot,  display,  display3d, 
fieldplot,fieldplot3d,  gradplot,  gradplot3d,  implicitplot,  implicitplot3d,  inequal,  listcontplot, 
listcontplot3d,  listdensityplot,  listplot,  listplot3d,  loglogplot,  logplot,  matrixplot,  odeplot,  pareto, 
pointplot,  pointplot3d,  polarplot,  polygonplot,  polygonplot3d,  polyhedra_supported, 
polyhedraplot,  replot,  rootlocus,  semilogplot,  setoptions,  setoptions3d,  spacecurve, 
sparsematrixplot,  sphereplot,  surf  data,  textplot,  textplot3d,  tubeplot ] 

[  arc,  arrow,  circle,  cone,  cuboid,  curve,  cutin,  cutout,  cylinder,  disk,  dodecahedron,  ellipse, 
ellipticArc,  hemisphere,  hexahedron,  homothety,  hyperbola,  icosahedron,  line,  octahedron, 
pieslice,  point,  polygon,  project,  rectangle,  reflect,  rotate,  scale,  semitorus,  sphere,  stellate, 
tetrahedron,  torus,  transform,  translate,  vrml] 
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%  THESIS 

%  PlotlT:  deltau  /  vc  vs .  v  /  vc  (or  Is  /  vc)  with  10%  TANKAGE  PENALTY 
%  "incremental  change  of  vehicle  velocity  vs.  impulse" 

%  LT  John  Jay  De  Beilis,  USN 
%  25  AUG  99 

%  Variables: 

%  x()  =  v  /  vc 

%  a ( )  =  mpl  /  mo 

%  y()  =  deltau  /  vc 

xl=linspace (0, 3, 300) 
x2=linspace (0,3,300) 
x3=linspace (0, 3, 300) 

al=0 . 35 
a2=0 . 45 
a3=0 . 55 

yl=xl . *log ((1.1  +  xl.A2)./(al  +  xl.A2  +  0.1)) 

y2=x2.*log( (1.1  +  x2.A2)./(a2  +  x2.A2  +  0.1)) 

y3=x3.*log( (1.1  +  x3.A2)./(a3  +  x3.A2  +  0.1)) 

figure (1) ,  plot (xl, yl, x2, y2, x3, y3) ,  grid 

axis( [0.0  2.0  0.0  0.5] ) 

ylabel ( 1 deltau  /  vc 1 ) , xlabel ( f v  /  vc’) 

gtext ( 'mpl/mo  =  0.351),  gtext ( 'mpl/mo  =  0.45*),  gtext ( 'mpl/mo  =  0.55’) 
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%  THESIS 

%  PlotlT :  deltau  /  vc  vs.  v  /  vc  (or  Is  /  vc)  with/without  10%  TANKAGE  PENALTY 
%  "incremental  change  of  vehicle  velocity  vs.  impulse" 

%  LT  John  Jay  De  Beilis,  USN 
%  25  AUG  99 

%  Variables: 

%  x ( )  =  v  /  vc 

%  a ( )  =  mpl  /  mo 

%  y()  =  deltau  /  vc 

xl=linspace (0, 3, 300) 
x2=linspace (0, 3, 300 )  , 
x3=linspace (0, 3, 300) 

a=0 .45 

yl=xl.*log( (1.1  +  xl . A2) . / (a  +  xl.A2  +  0.1)) 
y2=x2 . *log ( (1  +  x2.A2)./(a  +  x2.A2)) 

figured),  plot  (xl,  yl,x2,  y2) ,  grid 

axis( [0.0  2.0  0.0  0.5] ) 

ylabel (' deltau  /  vc' ) ,xlabel ( 'v  /  vc') 

gtext ( 'mpl/mo  =  0.45'),  gtext ( 'mpl/mo  =  0.45  with  10%  tankage  penalty') 
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